Self-avoiding knots 
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Scaling arguments are used to analyze the size of topologically constrained closed ring polymer 
with excluded volume. It is found that there exists a finite range of polymer thickness (excluded 
volume) in which self-avoidance is unimportant and polymer swelling compared to the Gaussian size 
is entirely due to the topological constraints. 

PACS numbers: 61.25.Hq, 82.45.-h, 82.20.Wt 
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The ability to knot is the most obvious and the 
least understood property of long strings, such as ropes, 
threads, spaghetti, flux lines, or polymers. Knots in DNA 
have been directly observed 0, El. The compelling, albeit 
circumstantial, proof of their importance can be seen in 
the fact that Nature had undertaken to design special cell 
machinery (called topological enzymes) spending energy 
to simplify DNA topology [||. Knots and other topolog- 
ical constraints are also of great importance in equilib- 
rium and dynamics properties of polymer materials, such 
as melts, gels, networks, etc. Despite this ample motiva- 
tion, current understanding of polymer topology is very 
limited. 

Apart from mathematical theory of knot invariants 
Q, the study of knots in polymers was dominated 
for a long time by computer simulations f^, [|. 

The incorporation of sophisticated topological invari- 
ants to compute entropy of the topologically constrained 
polymers remains a difficult challenge More re- 

cently, some insights into the statistical properties of 
knots were developed based on various scaling arguments 
0, [fL], [H| [13| PJllLAs in other cases in soft condensed 
matter physics |1 6|, |l7||, scaling is the theory of choice 
where more sophisticated approaches are too difficult. 

The simplest question of polymer physics is that about 
the size (say, root-mean-squared gyration radius) of a 
single macromolecule in a dilute solution. In general, 
this question is not understood for the topologically con- 
strained polymer ring. As a first step in this direction, it 
was recently shown ||12|] (confirming the earlier conjecture 
jt8|) that topological constraints in the knot effectively 
lead to self-avoidance on the large enough length scales. 
As a result, the size of a polymer ring with the topology 
of a trivial knot (an unknot) scales as R ~ tN v even for 
the polymer with negligible excluded volume. Here, N is 
the number of effective segments, I is the segment length, 
and v rs 0.588 ss 3/5 is the well known critical exponent 
of self-avoiding walks jl6| . This result has been also ex- 
tended for the non-trivial knots using the aspect ratio 
of maximally inflated tube, p, as a topological invariant 
IllOi. Overall, the result of the work |Ol reads 



Here, Nqo is the characteristic length of knotting; it 
comes from the observation [|| |tJ and mathematical proof 
p2j that the probability to have an unknot in a phantom 
chain (which freely crosses itself) goes like e~ N ^ N °°. Nu- 
merically, N Q0 ~ 300 [||, 0| . Thus, the theory predicts 
that R ~ A -3 / 5 , if N is large enough, independently of 
the knot type. This prediction was scrutinized using com- 
puter simulations Jli| [2(], |T) ; the most recent work H] 
claims complete agreement with the predicted scaling. 

In this paper, our goal is to study the interplay of 
topological constraints and excluded volume interactions 
- two factors both of which are inevitably present in any 
real polymer. We imagine a polymer ring of N effective 
segments of the length I and diameter d each. We would 
like to emphasize that end closing itself would not change 
the polymer size scaling, affecting only a numerical pref- 
actor. This means, if we imagine a ring with excluded 
volume, but phantom (for which all topological classes of 
conformations are available, like in the presence of topo- 
logical enzyme |2^]), then the usual prediction Jl(| holds: 
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IN 1 ' 2 when N < l/r 2 
m » r 2u-\ when TV > l/r 2 
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Thus, our goal if to generalize together the Eqs. ([!]) and 
(^). Our results are summarized by the Fig. [l], in which 
the diagram of scaling regimes is presented in terms of 
variables N and r = d/l. The paper is organized around 
the discussion of this diagram. 

As an important input of our analysis, we shall rely 
on the following insight gained by computer experiments 
pi] , psf regarding the probability of knotting for the ex- 
cluded volume polymer. As one might have guessed, 
excluded volume enhances the probability of a trivial 
knot, somewhat suppressing all non-trivial knots. What 
is much more difficult to guess is that the probability 
of a trivial knot remains exponential in N, with charac- 
teristic length very sharply depending on the excluded 
volume parameter. According to 0, |25|| , 
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In fact, the authors of the work [g5| found that according 
to their numerical data, the value of v is slightly below 
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FIG. 1: Diagram of scaling regimes of a polymer ring. N is 
the number of effective segments in the chain, r is the dimen- 
sionless thickness of the polymer (the ratio of the diameter to 
the segment length) , p is the knot invariant - diameter to the 
length ratio of maximally inflated tube [ju^]. The scaling of 
polymer size for every regime is indicated in the figure using 
the approximate value of the critical exponent v w 3/5. The 
cross-over boundaries are as follows: the G — S line between 
Gaussian and self-avoidance regimes r 2 ~ p/N; the G — T 
line between Gaussian and topological regimes N ~ pNo(r) 
(see Eq. ([])); and the T — S line between topological and 
self-avoidance regimes r = r* « 0.03. 



unity, v « 0.85. They also claimed that the coefficient 
A me depend on the model. Neither of these details is 
important for our purposes here, we will only use the 
fact that the characteristic length of knotting, No(r) is 
sharply increasing with r. 

We are now in a position to start discussing the scaling 
diagram, Fig. |l|. Let us begin with a trivial knot, or an 
unknot, in which case p (in the scaling sense) is of order 
unity (p = 271"). 

First of all, there is obviously a regime when neither 
excluded volume nor topological constraints are impor- 
tant. This range is labeled as Gaussian in Fig. |], because 
chain size in this regime scales as that of a Gaussian poly- 
mer, R ~ N 1 / 2 . This happens if the conditions N < 1/r 2 
(||) and N < N (r) are met simultaneously. In par- 
ticular, the latter condition, according to the formula (||) 
ensures that even for the phantom chain the part of con- 
formation space corresponding to all non-trivial knots is 
very small; therefore, set of conformations of phantom 
chain is practically the same as that of a trivial knot. 

When we cross over the G — S boundary N ~ 1/r 2 , 
excluded volume effect takes over, but topological con- 
straints may still remain unimportant if N < No(r). In 
this regime, labeled as dominated by self-avoidance in 
Fig. |l[ chain size is given by R ~ tN v r 2v ~ x . 

Another regime is obtained from Gaussian if we cross 
over the N = No(r) line G — T, remaining at N < 1/r 2 . 
In this case, we should imagine the chain as consist- 



ing of N/No(r) blobs, of No(r) monomers each. Inside 
each blob, as it length is smaller than the characteris- 
tic knotting length, the topological constraints are in- 
significant, as well as the excluded volume; therefore, 
the blob size is about i?biob — (.N^r) 1 ^ 2 . On the other 
hand, the chain of blobs is topologically constrained, and, 
therefore, according to fl2fl , we get for the chain size 
R ~ i? blob {N/N {r)) v ~ £N v N (r)-" +1 / 2 . 

At this stage, it is useful to check the self-consistency 
of our result. As we know 10, [TJ|], the trivial knot size 
results from the balancing of two entropic factors. On the 
one hand, chain has entropic elasticity ~ R 2 /Ni 2 which 
resists swelling. On the other hand, topologically con- 
strained blobs resist interpenetration. At the equilibrium 
R, both entropic contributions are of the same order. We 
must now check that the second virial term correspond- 
ing to the excluded volume interactions is smaller than 
either of the above mentioned terms. Given that the sec- 
ond virial coefficient of elongated (rod-like) segments is 
about £ 2 d, the self-consistency condition reads 
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Note, that this condition is equivalent to R > i?Fiory, 
where i?Fiory — IN^^r 1 ^ is the Flory estimate of poly- 
mer size without topological constraints [ ^6[ . Thus, this 
condition means simply that excluded volume is unim- 
portant if the chain swells more for topological reasons 
than it would due to the excluded volume. Simple alge- 
bra indicates that this condition is met if No(r) < 1/r 2 , 
which is equivalent to r < r* . Numerically, according to 
eq. (§, r* » 0.03. 

Finally, let us see what happens if No(r) > 1/r 2 
and N > No(r). Diagram Fig. [I] indicates that in 
this case we cross over the T — S boundary and en- 
ter again the same self- avoidance regime which we dis- 
cussed previously n terms of crossing over the G — S 
line. Let us show why this is indeed the same regime. 
We should imgine once again that our chain consists 
of N/No(r) blobs, with No(r) monomers in each blob. 
In the present case, since No(r) > 1/r 2 , excluded vol- 
ume is important for each blob, so that the blob size is 
about i?biob — INoiryr 2 "' 1 . Since the chain of blobs 
is topologically constrained, overall chain size scales as 
R ~ Ruob {N/N {r)) u ~ ^Vr 2 " -1 , which does not in- 
volve the topological quantity of No (r) and coincides with 
our result for the self-avoidance dominated regime. 

This completes the analysis for the case of an unknot 
(trivial knot). 

Let us now consider some non-trivial knot, character- 
ized by the length-to-diameter ratio of maximally inflated 
tube, also called ideal knot representation. The specific 
values of p for many knots were computed by Pieranskii 
et al [^6| . Following the ideology of 10 [l2| , we imagine 



the polymer as confined in the maximally inflated tube 
such that it is unknotted inside the tube. This corre- 
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sponds to viewing the polymer as consisting of p roughly 
spherical blobs, with N/p monomers in each blob. Fur- 
thermore, since the tube is maximally inflated, the blobs 
fill completely (in the scaling sense) the entire volume of 
the coil, which means R 3 ~ pi?j^ lob , or 



R ~ p 1//3 i?blob • 



(5) 



Therefore, all we have to do is to address all possible 
regimes of the chain inside the blob. 

The simplest case is when neither excluded volume nor 
topology are of importance inside the blob. This happens 
when both N/p < l/r 2 , which is below the G — S line in 
Fig. [I], and N/p < N a (r), which is above or to the left 
of the G — T line. In this regime every blob is Gaussian, 
yielding i?biob — ^{N/p) 1 / 2 and, according to formula 



i?~ £N 1/2 p- 1/& ,when N < p/r 2 

and N < pN (r) 



This is Gaussian regime of Fig. |[ 

When N/p exceeds l/r 2 , we cross over the G — S line 
into the regime dominated by self-avoidance. In this case, 
the blob size is about i? Wo b ^ l(N / p) v r^ 1 . Therefore, 
formula (||) yields 

R ~ m u p- y+1/z r 2y - 1 , when TV > p/r 2 

and r > r* . (7) 

Similar to the case of the trivial knot, this regime, which 
we call self-avoidance dominated, is valid all the way from 
the G — S boundary N ~ p/r 2 to the T — S boundary 
r = r* . 

Finally, if r is small, but the blob length N/p exceeds 
No(r), then blobs become topologically constrained. This 
means, it becomes important that the chain inside the 
tube is unknotted, and so must be every blob. In this 
regime, blob size scales as i?biob — £(N/pYN (r)~ u+1 / 2 . 
Plugging this into formula (||) leads to 

i?~ lN"p- u+1/3 N (r)-' /+1/2 ,when N > pN (r) 

and r < r* . (8) 

This is topology dominated regime. 

It is instructive to summarize the results by looking at 
the ratio of the given knot size R to the size of a phantom 
ring polymer (Q), with the same length N, segment length 
L and excluded volume r: 



-Rknot(p) 
-^phantom 



(9) 



This quantity is useful, because -Rphantom can be also un- 
derstood as the average size over all knot types, and such 
quantity can be easily extracted from computer simula- 
tions data [jl9L po). Various regimes for this quantity are 
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FIG. 2: Scaling regimes for the quantity £ (|9|). As in Fig. 
[j], the results are shown using the approximate value of the 
critical exponent v ~ 3/5. Regime A: both topologically con- 



strained and phantom chains are Gaussian, £ ~ p 



-1/6 



< 1. 



(6) Regime B: phantom polymer swells due to excluded volume, 



while knot remains Gaussian: £ 



1/6 



(Nr 2 



-1/10 



< 1. 



Regime C: both phantom polymer and the knot swell due to 
the excluded volume, £ ~ p _4//15 < 1. Regime D: knot swells 
because of topology, while phantom ring swells due to the ex- 
Regime E: knot is 



eluded volume, £ ~ p~ i/15 (N (r)r 2 ) 



-1/10 



swollen due to the topological constraint, while phantom ring 
is Gaussian, £ ~ p" 4/15 (N/N (r)) 1/10 . In the shaded area, 
which includes parts of the regimes D and E, the topologically 
constrained polymer is more expanded than the phantom (av- 
eraged of all knots) . 



outlined in the Fig. E, where the cross-over lines of the 
scaling diagram Fig. I] are superimposed with the line 
r ~ 1/yN separating Gaussian and swollen regimes for 
the phantom polymer. 

Somewhat counterintuitive result is that there is a pa- 
rameter range (shaded in Fig. ||) in which topologically 
constrained knot is more expanded than phantom aver- 
age over all knots. This happens because for the phan- 
tom chain (with annealed topological disorder) in this 
range more complex knots (with p larger than the given 
value) are prevailing, they are more compact, and they 
dominate the average -Rphantom- This result is at least 
qualitatively in agreement with the simulation data p9| . 

Metzler et al jl3|, |l4| recently reexamined the possibil- 
ity of "knot segregation," which was briefly discussed at 
the end of |[o| and later analyzed numerically (^7). The 
idea is that the entropy of a knotted polymer may be 
maximaized by the conformations in which knot is tight- 
ened within a short part of the chain, while the rest of the 
polymer fluctuates as a free unknotted loop. While the 
tendency to tighten the knot is quite strong for charged 
polymers |]l5| , for neutral system it is governed by sub- 
linear in TV terms of free energy and, therefore, is very 
sensitive to the details of structure and interactions. In 
this paper, we assumed that topological constraints are 
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delocalized along the polymer chain backbone. In the 
context of our work, this is justified for the following rea- 
sons. First, this issue is altogether irrelevant for the im- 
portant case of the trivial knot, which is the first case con- 
sidered in this paper. Second, for the non-trivial knots, 
it is necessary to understand the delocalized states in or- 
der to be able to push the question of knot segregation 



beyond the slip- link model of 14 



It should also be noted that in this paper we always 
assumed knots simple enough such that polymer chain is 
far from fully stretched state in the tube. The regime of 
such strong knotting requires special analysis |2§1| . 



Tt> conclude, we presented the scaling theory determm- 
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ing the size ot a polymer ring with quenched topology ot 
a certain knot. 
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